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Abstract 

We show that in a large class of models based on anomalous U(l) symme- 
try which addresses the fermion mass hierarchy problem, leptonic flavor changing 
processes are induced that are in the experimentally interesting range. The fla- 
vor violation occurs through the renormalization group evolution of the soft SUSY 
breaking parameters between the string scale and the ^(1)^ breaking scale. We 
derive general expressions for the evolution of these parameters in the presence of 
higher dimensional operators. Several sources for the flavor violation are identi- 
fied: flavor-dependent contributions to the soft masses from the U(1)a gaugino, 
scalar mass corrections proportional to the trace of U(1)a charge, non-proportional 
A-terms from vertex corrections, and the U(1)a D— term. Quantitative estimates 
for the decays fi — > e~f and t — ► /ry are presented in supergravity models which 
accommodate the relic abundance of neutralino dark matter. 
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1 Introduction 



The Standard Model, while highly successful in explaining all experimental data, does not 
provide an explanation for the observed hierarchy in the masses and mixings of quarks 
and leptons. Extended symmetries are often speculated to address this problem. Family- 
dependent U(l) symmetry is a widely studied extension. An attractive scenario is the 
Froggatt-Nielsen scheme pQ. In this scenario all the Yukawa couplings are assumed to 
be of order one, but the ones which generate the light fermion masses arise only as 
nonrenormalizable operators suppressed by powers of a small parameter e = (S)/M, 
where (S) is the flavor symmetry breaking order parameter and M is a more fundamental 
mass scale. With the flavor U(l) charges of fermions differing only by order one, large 
hierarchy factors, such as m u /m t ~ 10~ 6 , are explained. 

A natural origin for the flavor U(l) symmetry is the anomalous U(1)a gauge symmetry 
of string theory [2] . The small expansion parameter e arises in a natural way in anomalous 
£7(1) models through the Fayet-Iliopoulos term induced by the gravitational anomaly [HJ. 
Such models have been extensively studied in the literature for understanding the fermion 
mass hierarchy puzzle [HE]. The purpose of this paper is to address flavor changing neutral 
currents in this class of models. 

In the presence of low energy supersymmetry (SUSY), a family-dependent anomalous 
U(1)a symmetry will induce flavor changing processes, even when there is no flavor vio- 
lation in the soft SUSY breaking parameters at the fundamental scale - taken to be the 
string scale. Such violations will be generated through the renormalization group evolu- 
tion (RGE) of the SUSY breaking parameters between M string and the U(1)a breaking 
scale (S). We derive general expressions for the evolution of these parameters in the 
presence of higher dimensional operators. Our results can be applied to a wide class 
of Froggatt-Nielsen models. We have found several sources of flavor violation. In the 
momentum interval (S) < [A < M string , the U(1)a gaugino is active and will contribute 
differently to the soft masses of different families. Because TtU(1)a is not zero in anoma- 
lous U(l) models, there are nonuniversal RGE contributions to the soft scalar masses 
arising from the D-term proportional to the respective flavor charges. Furthermore, the 
trilinear A-terms will receive vertex corrections from the U(1)a gaugino that are not 
proportional to the respective Yukawa couplings. In addition to these RGE effects, upon 
symmetry breaking, the D-term associated with the U(1)a will also induce nonuniversal 
masses for the sfermions |B] [!]• 

In this paper we investigate the combined effects of nonuniversality for lepton flavor 
violating (LFV) decays \i — > and r — > ^7 in a class of anomalous U{1)a models. 
Quantitative predictions for the branching ratios are presented in two specific models 
of fermion mass hierarchy. We find that the branching ratio for /1 — ► ej is around the 
current experimental limit, while r —>■ fij may be accessible in the future. In our analysis 
we also include the right-handed neutrino-induced LFV effects, which have been widely 
studied in the literature jSj-[H]. These effects turn out to be significant in some but not 
all cases that we study. In fact, within our framework, if leptogenesis is assumed to be the 
source of the observed cosmological baryon asymmetry , it turns out that the right- 
handed neutrino induced LFV effects are negligible. Related effects from SUSY GUT 
thresholds for LFV have been studied in Ref. |13j . The flavor violating effects are more 
prominent in the leptonic sector compared to the quark sector since quark flavor violation 
is diluted somewhat by the gluino focusing effects which arise during the evolution of the 
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SUSY breaking parameters below the U(1)a breaking scale. This is especially so when 
the SUSY parameter are chosen, as we do, such that the lightest neutralino is the dark 
matter with an acceptable cosmological abundance. 

The structure of the paper is as follows. In Section 2 we introduce our models. In 
2.1 we describe our models of fermion mass hierarchy, in 2.2 we present the details of 
anomalous U(l) models. In 2.3 we present our fermion mass fits for two models, Model 
1 and Model 2. In 2.4 we discuss baryogenesis through leptogenesis in these models. In 
Section 3 we derive generalized expressions for the evolution of the soft SUSY breaking 
parameters appropriate for a class of Froggatt-Nielsen models. Section 4 is devoted to the 
numerical analysis of the branching ratios for the processes \x — > ej and r — > /ry. In 4.1 
we outline the qualitative features of flavor violation arising from various sources. 4.2 has 
our numerical fits. Figures 3 and 4 show our results for the branching ratios for fx — > 
in the two models of fermion masses. In 4.3 we address the D-term splitting problem. 
Our conclusions are given in Section 5. 



2 Fermion Mass Matrices and Anomalous U(1)a Fla- 
vor Symmetry 

In this section we present two specific models of fermion mass hierarchy based on a 
flavor-dependent anomalous U(1)a symmetry using the Froggatt-Nielsen mechanism pQ. 
The U(1)a symmetry is broken by an MSSM singlet flavon filed S slightly below the 
string scale (M st ). This provides the small expansion parameter e = (S)/M st needed 
for explaining the fermion mass hierarchy. We also present explicit models of anomalous 
U{1)a symmetry. 



2.1 Fermion Mass Hierarchy 

A general form of the superpotential which can explain the fermion masses and mixing 
hierarchy through the Froggatt-Nielsen mechanism has the form 



Vii„ f s \ n % , yt<„ ( s ^ " 
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n-ji ' \M st J <,! 1 J u \M st 
+ M Rij ^(—) V +nH u H d , (1) 
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where i,j = (1, 2, 3) are family indices, nfj, nfj, nf^, and are positive integers fixed 
by the choice of U(1)a charge assignment, yfj etc. are Yukawa coupling coefficients which 
are all taken to be of order one. Here Mr is the right-handed (u£) neutrino mass scale. 
We use the standard notation for the MSSM fields. 

The soft supersymmetry breaking terms which will induce LFV have the form given 

by 
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- Mi /2 X i X i + -M XF X F X F + B^H u H d + h. c. 
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+ E^Ui+^W+W+W- ( 2 ) 

Here a tilde stands for the scalar components of the matter superfields, and Aj and Xp are 
the MSSM gauginos and the flavor U(1)a gaugino. (My 2 , M\ F ) and ((mj) a b, fh 2 s ) are the 

gaugino and scalar soft masses respectively. f a stands for the MSSM sfermions including 
the right-handed sneutrinos. Note that the generalized A-terms in Eq. (J2J) has the same 
structure as the corresponding superpotential terms of Eq. (JTJ). 

We assign flavor U(1)a charges to the MSSM fields such that the observed fermion 
mass and mixing hierarchies are obtained with all Yukawa couplings being order one. As 
we will show explicitly, the expansion parameter e = {S)/M st is naturally of order 0.2 in 
anomalous U(l) models. We use the idea of "lopsided" mass matrices for generating large 
neutrino mixings fUE3, while maintaining small quark mixings. 

We analyze two specific textures for the fermion mass matrices, Model 1 and Model 
2. Model 2 is the texture advocated in Ref. [5], and Model 1 is a slight variant of it. The 
textures for the two models are 
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Here M u , M d and M e are the up-quark, down-quark, and the charged lepton mass ma- 
trices (written in the basis uM u u c , etc.). M UD is the Dirac neutrino mass matrix, and 
M v c is the right-handed neutrino Majorana mass matrix. The light neutrino mass matrix 
flight j g d er i ve( } f r0 m the seesaw mechanism. We have not exhibited order one coefficients 
in the matrix elements of Eq. (J3*|). The expansion parameter is e ~ 0.2. The exponent 
p appearing in the overall factor e p multiplying M d and M e is assumed to take values 
0, 1 or 2 corresponding to large (~ 20), moderate (~ 10), and small (~ 5) values of 
tan/? (= (H u ) / (H d )) respectively. 

The parameter a is allowed to take two values, and 1, corresponding to Model 1 
(a = 0) and Model 2 (a = 1). The two models differ only in the masses and mixings of 
the first family. Both models give excellent fits to the fermion masses and mixings includ- 
ing neutrino oscillation parameters. Their predictions for LFV are however noticeably 
different, which we analyze in Section 4. 



2.2 Anomalous U{l) A Models 

The mass matrix textures of Eq. Q can be obtained from an anomalous U(l) gauge 
symmetry that is flavor dependent. In Table H we present a set of such £7(1) charges. 
These charges are integer multiples of the charge of the flavon field S which is taken to 
be —1. Note that the texture alone does not fix the exponent s appearing in M UD of Eq. 
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H u , Hd, S 


0, 0, -1 


(h,h,q s ) 



Table 1: The flavor U(1)a charge assignments for the MSSM fields and the flavon field S in the 
normalization of q s = —1. Here a is (1) for Model 1 (Model 2). In the third column we list the generic 
notation for the charges used in the RGE analysis. 



We use the Green-Schwarz (GS) mechanism |2j for anomaly cancellation associated 
with U(1)a gauge symmetry. These conditions are given by 



A\ _ A 2 _ A 3 _ Aa _ A gravity 



(4) 
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where M and A A are the coefficients U{\)\ x U(1) A , SU (2)1 x U(l) A , SU(3) 2 C x U(1) A 
and U(1) 3 A gauge anomalies respectively, ki and k A are the Kac-Moody levels. A gravity 
is the mixed gravitational anomaly coefficient which is given by the trace of the U(l) a 
charges over all fields. With the non-Abelian levels h% — ks — 1, which is the simplest 
possibility, from Table [T] and Eq. (jlj) one finds 



A, 



19 - 3a + 3s 
2 : 



. 19 - 3a + 3p . . 

^3 = ^ ■ (5) 

This implies that p = s. Furthermore, 

A x = ^(19-3a + 3p), (6) 
o 

which fixes the level k\ to be 5/3. 

With p = s the charges given in Table d are compatible with SU(5) unification. In 
string theory gauge coupling unification can occur without a simple covering group. The 
string unification condition is ^H] 

hgf = k 2 g 2 2 = hgl = k A g 2 F . (7) 
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Our result k\ = 5/3 is what is needed for consistency of the observed unification of gauge 
couplings in the MSSM. The discrepancy in the unification scale derived from low energy 
data versus perturbative string theory evaluation can be reconciled in the context of M- 
theory by making use of the radius of the eleventh dimension [Fj\. We assume such a 
scenario. 

From now on we shall assume the SU(5) normalization for g\. If one assumes that 
the field content of the model is just the one listed in Tabled the gravitational anomaly 
-^gravity would not satisfy the GS condition. One simple solution is the introduction of 
additional MSSM singlet (hidden sector) fields X^. Then Eq. (J2J) leads to the following 
result: 

A gramt y = 5 (13 - 2a + 3p) + £ gjf = 6(19 - 3a + 3p) , (8) 

k 

where q x are the U{1)a charges of the extra fields X k . From this, one gets J2k1k = 
(49 — 8a + 3p) . We assume for simplicity that all the X k fields have the same flavor 
charge equal to 1. The number n x of X k fields is then fixed to be 

n x = 49 - 8a + 3p . (9) 

We are now in a position to determine the level as well as the U{1)a gauge coupling 
gp at the unification scale. We renormalize the U(1)a charges by a factor \q s \ so that the 
charge of the flavon field is now — |<7 S |. \q s \ is determined by demanding gp = g\ at the 
unification scale. Eq. (0J) and the number n x in Eq. Q then fix \q s \ to be 



19 - 3n -i Hp 



x 



V 10(4-a) 3 + 5[(l+p) 3 + 2p 3 ] +n- 

For p = (0, 1, 2) one has \q s \ = (0.165, 0.172, 0.166) for Model 1 (a = 0) and \q s \ = 
(0.225, 0.228, 0.203) for Model 2 (a = 1) . 

The Fayet-Iliopoulos term for the anomalous U(1)a, generated through the gravita- 
tional anomaly, is given by [3] 



, _ 9 St M St I ] | , | | , 

^ — 2 Ws]-™- gravity , - 11 > 



where g s t is the unified gauge coupling at the string scale. By minimizing the potential 

v = hM I e _ |5|2 + E q f lIal 2 + E q x lx A 2 (12) 

8 \\Qs\ a k J 



in the unbroken supersymmetric limit we obtain 

(13) 



(S) 1 9st A gravity 



M st V 192tt 2 

The numerical values of e derived from Eq. (|13p for different p and a are listed in Table 121 
by making use of Eq. ()10|) . This is the small expansion parameter appearing in the mass 
matrices of Eq. (jHJ). Here we took g^j^K ~ 1/24. 
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e 


p = 


p — 1 


p = 2 


a = 


0.177 


0.191 


0.204 


a = 1 


0.163 


0.177 


0.191 



Table 2: Numerical values for the small expansion parameter e corresponding to different fermion mass 
hierarchy structure. 



The mass of the flavor U(1)a gauge boson is found to be 

M F = lqsl9 f } . (14) 
v2 

Between the string scale M st and Mp the flavor gaugino contributes to flavor violating 
processes. This mass can now be determined: 

V88.3' 84.6' 86.3/ V ' K ' 

in the case of Model 1 and 

M st M st M st 



in the case of Model 2. 



68.7' 66.8' 72 



^f) for p= (0,1, 2), (16) 



2.3 Fermion Mass Fits 

Here we present numerical fits to the fermion masses and mixings for Model 1 and Model 
2. These fits will be used in our quantitative analysis of lepton flavor violation. 

As input at low energy we choose the following values for the running quark masses 

CHI 

m M (lGeV) = 5.11 MeV, m c (m c ) = 1.27 GeV, m t (m t ) = 167 GeV, 
m d (l GeV) = 8.9 MeV, m s (l GeV) = 130 MeV, m b (m b ) = 4.25 GeV. (17) 

The CKM mixing matrix elements are chosen to be |V^ S | = 0.222, \V u b\ = 0.0035, \V c b\ = 
0.04 and rj = 0.33 (the Wolfenstein parameter of CP- violation) . Using two-loop QED 
and QCD renormalization group equations we obtain these running parameters at the 
SUSY breaking scale, M SUS y = 500 GeV, with a s (M z ) = 0.118, to be 

r ^m f (M SUS Y) , 1C , 



m f {m f ) 
where 

(r t , r b , r T , r UjC , r d ^ r e , M ) = (0.943, 0.605, 0.991, 0.395, 0.398, 0.989) . (19) 
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Using two-loop SUSY RGE evaluation above M S usy we obtain the Yukawa couplings at 
the U(l) A breaking scale (~ 10 15 GeV) to be 

(Y u , Y c , Y t ) = (5.135 x HT 6 , 1.426 x 10~ 3 , 0.538) , 

(Y d , Y s , Y b ) = (3.459 x 10~ 5 , 5.052 x 10" 4 , 2.768 x 10~ 2 ) , 

(Y e , Y^ Y T ) = (1.024 x 10~ 5 , 2.118 x 10~ 3 , 3.572 x 10~ 2 ) , 

(Y U1 , Y„ 2 , Y U[i ) = (3.515 x 10~ 4 , 8.419 x 10~ 4 , 1.131 x 10~ 2 ) , (20) 

for tan (3 — 5, 

(Y u , Y c , Y t ) = (4.999 x HT 6 , 1.389 x HT 3 , 0.518) , 

(Y d , Y s , Y b ) = (6.844 x 10~ 5 , 9.997 x 10~ 4 , 5.470 x 10~ 2 ) , 

(Y e , Y^ Y T ) = (2.027 x 10~ 5 , 4.192 x 10~ 3 , 7.094 x 10~ 2 ) , 

(Y Ul , Y U2 , Y V3 ) = (1.708 x 10~ 3 , 4.105 x 10~ 2 , 5.519 x 10~ 2 ) , (21) 

for tan (3 = 10, and 

(Y u , Y c , Y t ) = (4.996 x 10~ 6 , 1.387 x 10~ 3 , 0.518) , 

(Y d , Y s , Y h ) = (l.40 x 10~ 4 , 2.045 x 10~ 4 , 0.113) , 

(Y e , Y^, Y T ) = (4.132 x 10" 5 , 8.545 x 10~ 3 , 0.147) , 

(Y V1 , Y U2 , Y„ 3 ) = (8.551 x 10" 3 , 2.059 x 10" 2 , 0.278) , (22) 

for tan/? = 20. \V u b\, \V c b\, \V t d\ and \V ts \ are multiplicatively renormalized by an RGE 
factor of 0.9 in going from the low energy scale to the U(1)a breaking scale. 

We have determined the Dirac neutrino Yukawa couplings as follows. First we note 
that the anomaly cancellation conditions in Eq. (j3J) implies p = s, which means that the 
Dirac neutrino Yukawa couplings are fixed to be of the same order as the charged lepton 
Yukawa couplings. Now, if one takes the right-handed Majorana neutrino mass matrix 
to be proportional to the transpose of the Dirac neutrino Yukawa coupling matrix for 
simplicity, M u c = Yj Mr e p , then the light neutrino mass matrix is given by 

M ij gM = YvM - }Y T v 2 sm 2 p = yfj^ll . (23) 

M R eP 

This simplified choice is certainly consistent with the fermion mass structures we have 
chosen in Eqs. (JHJ). We adopt this choice in our analysis. Y v is determined from a fit to 
the light neutrino oscillation parameters with Mr = 10 14 GeV. This fit corresponds to 
m Ve = 2.7 x 10~ 3 eV, m V( , = 6.4 x 10~ 3 eV and m Ur = 8.6 x 10 -2 eV and the leptonic 
mixing matrix given by 

/ 0.848 -0.526 -0.0409 \ 
V MNS = 0.349 0.619 -0.72 . (24) 
V-0.4 -0.59 -0.7013/ 
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We also consider a scenario where the Dirac neutrino Yukawa couplings are maximized 
by choosing M« = 4 x 10 14 GeV. In this case we have 



(1.406 x HT 3 , 3.368 x HT 3 , 0.453 x 10~ 2 ) for tan/3 = 5, 

(6.843 x 10" 3 , 1.645 x HT 2 , 0.222) for tan/3 = 10, 

(3.514 x 10" 2 , 8.464 x 10~ 2 , 1.237) for tan/3 = 20. (25) 



We now present our fits to the observables of Eqs. pU jl - lJ^Sjl consistent with the texture 
of Eq. ©. This cannot be done uniquely since the right-handed rotation matrices are 
unknown from low energy data, so we make a specific choice. In our lepton flavor violation 
analysis we shall make use of this specific fit. One should bear in mind that there are 
uncertain coefficients of order one in the Yukawa matrices of our fit, which can lead to an 
order of magnitude uncertainty in the branching ratios for LFV processes. 

We introduce the following notation: 



n 



f n 1 . 



(26) 



In Model 1, a good fit to the Yukawa couplings matrices is found to be 



3/33 



/ 3.91 e 8 
0.226 e 6 
V 0.375 e 4 



0.226 e 6 
1.91 e 4 
0.499 e 2 




(1.56 + 0.115t) e 5 
-2.89 e 3 
C— 0.878 + 0.88 x 10 3 ? 



(0.909 + 0.054z)e 4 (0.658 + 0.131i) e 4 
1.02 e 2 1.22 e 2 

le 0.412 + 0.11 x 10~ 6 i 1 + 0.73x10^ 



ye 



f 1.89e 5 1.57e 3 0.812e^ 

0.487e 4 2.14e 2 0.316 
V l.lOe 4 1.52e 2 1 



Here 



( 1.51 e 2 
-0.358 e 
-0.438 e 



-0.358 e 
0.339 
0.485 




(27) 



2/33 
2/33 
2/3 £ 3 
2/33 



0.539, 0.523, 0.519) 
0.650, 0.257, 0.106) 
0.840, 0.333, 0.139) 
0.225, 0.219, 0.221) 



(2? 



for (p = 2, 1, 0) which we shall associate with tan/3 = (5, 10, 20). Here we have taken 
e = 0.2. For simplicity we assumed the leptonic Yukawa couplings to be all real. 
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In Model 2 we have the following fit for the Yukawa coupling matrices: 



Jo 



2/33 



f 0.876 

1.30 e 5 
V 0.499 e E 



1.30 e 5 
2.59 e 4 
0.993 e 2 



0.499 e 3 ' 
0.993 e 2 
1 



I (3.01 + 0.13i)e 4 

1.79 e 3 
V (1.00 + 0.33 x nT 3 i)e 



(2.66 + 0.13i)e 3 
2.26 e 2 
0.987 + 0.582 x 10 



-5 • 



(1.21 + 0.13z)e 3 

1.42 e 2 
1 + 0.21 x 10 5 z 



/ 1.19 e 4 
ye = 2/33 eP 0.892 e 3 
V 1.36 e 3 



1.68 e 3 0.579 e\ 
2.18e 2 0.350 , 
1.45 e 2 1 / 




Here 



1.53 e 2 -0.329 e 
Y v = y% 3 e p | -0.329 e 0.293 0.449 j. (29) 
-0.406 e 0.449 



y% 3 = (0.535, 0.52, 0.515), 
y ^ = (0.650, 0.257, 0.107), 
y e 33 = (0.840, 0.333, 0.139), 

yl z = ( 0.233, 0.228, 0.229) , (30) 



for three different values of p = (2, 1, 0) identified with tan/3 = (5, 10, 20). 



2.4 Leptogenesis Constraints 

In this subsection we show that baryogenesis through leptogenesis jT2, can occur naturally 
in our models of fermion mass hierarchy. This requirement does put significant restrictions 
on the parameter p, which is related to the value of tan/3. We find that an acceptable 
baryon asymmetry would require p — 2, 1, corresponding to small to moderate values of 
tan p. In this case, lepton flavor violation resulting from the right-handed neutrino Dirac 
Yukawa couplings tend to be suppressed. 

From the neutrino Majorana mass matrix of Eq. (3), we find that the heavy Majorana 
masses of denoted by M 2j 3 are of order M R , while that of v\ (denoted by Mj) is of 
order e 2 Mr. Lepton asymmetry is generated in the out of equilibrium decay of v\. The 
induced asymmetry parameter is given by 




where v u ~ 174 sin/? GeV, and the function f(x) ~ —3/y/x. Let us focus on the asym- 
metry induced through the exchange of v%. Denote the mass ratio Mi/M 2 = ae 2 , where 
a is a coefficient of order one. We also have from Eq. (3) {Ml M UD ) U = be 2s+2 v^, 
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(Ml D M UD )i2 = ce 2s+1 vl, where b, c are also order one coefficients. e\ can then be esti- 
mated to be 

e i-|^+ 2 sin0, (32) 
8nb 

where sin0 is an order one phase parameter. The leptonic asymmetry parameter Yi is 
obtained from e\ via the relation Yl = de\/g*, where d is the dilution factor and g* is the 
number of degrees of freedom in thermal equilibrium (g* ~ 200 in the supersymmetric 
Standard Model). The dilution factor d depends on the ratio k of the decay rate versus 
the Hubble expansion rate |19j : 

Mr,, ( A/ft /If.. \ 

(33) 

We now demand that the light neutrino mass m^, which can be written from the last 
expression of Eq. Q as m U3 = he 2s v 2 /M 2 , should equal 0.05 eV, to be consistent with the 
atmospheric neutrino oscillation data. This then determines k to be k ~ [b/(ah)](33.3) . 
For this range of k (assuming that b/ (ah) is not smaller than about 0.3) the dilution factor 
can be written approximately as 

d^ , ri ;* R , (10<fc< 10 6 ). (34) 



M Pl (Mj D M v 
1.66^(8^2)71^1 ' 



Using Y B ~ —Y L /2, we find 



A;[ln(fc)] - 6 ' 



_ 3.6 x 1Q-<W 

K33.3y)]«-« ' 1 ' 



where x = (a 2 c 2 h sin <p)/b 2 and y = b/(ah) are both expected to be of order one. For the 
case of s = 2 , we choose x = 0.05, y — 3 to obtain Y B = 1 x 10~ 10 , which is near the 
central value of the observed baryon asymmetry. Note that this choice of x, y is quite 
natural and consistent with the flavor structure we have adopted. For the case of s — 1 , 
we can fit Y B = 1 x 10~ 10 by choosing x = 1.7 x 10~ 3 , y = 3 . Considering that x is 
a nontrivial combination of order one parameters, this choice of x cannot be considered 
unnatural. On the other hand, for the case of s = , a good fit to Y B requires, for 
example, x = 7 x 10~ 5 , y = 3. Such a low value of x does not go well with the hierarchical 
structure dictated by the anomalous U{1)a symmetry. 

We conclude that baryogenesis via leptogenesis works in a simple and natural way 
with the mass matrix textures suggested in Eq. (3), but only for low to medium values 
of tan (3 , corresponding to s = p = 1, 2 . 



3 Generalized RGE Analysis of Soft SUSY Breaking 
Parameters 

In this section we give a general RGE analysis of the soft SUSY breaking parameters 
including higher dimensional operators as shown in Eqs. (0) and (0). This includes the 
effects of the flavor U(1)a gaugino sector. Our analysis of this section should apply to a 
large class of Froggatt-Nielsen models. 
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The one-loop /3-functions for the soft scalar masses of the sleptons are found to be 

MSSM 1 



'J 



+2 (Y^mlY v + m 2 H Y v] Y v + A^A") 

+2qfg 2 F 5 l3 (a-Aqf(M XF ) 2 )} , (36) 

P (ml).. = /3(m^ /S5A %^ 2 g^.^_ 4g e (MAF) 2^ (3?) 

+4 (Vm 2 !" 4 + m 2 H Y v Y u] + A" A" 4 ) 

+2gr^(a-4gr(M AF ) 2 )}. (38) 

Similarly the /^-functions for the squark soft masses are given by 

(3 (ml),, = /3(my^ SA %^2g^^a-4g«(M AF ) 2 ), (39) 

P(m 2 u ) t . = (3 (m 2 u )* ISSM + ^2g^(a-4gr(M AF ) 2 ), (40) 



(3 (m 2 ) = (3 (m 2 )f SM + ^qfg^ (- - 4gf (M Af ) 2 ) . (41) 



Here a is defined as 



a = 3 Tr (2g Q m 2 3 + q u m 2 u + /m 2 ,) + Tr (2q L m 2 L + q e m\ + q v m 2 v 

+Qsm 2 s +Y / Qkrh 2 Xk , (42) 

A; 

where rhx k is the soft mass of the extra particles and the trace is over family space. 
Here j3 (w 2 )^ stands for the MSSM /3-function without the v c or the flavor U(1)a 
contributions |2Uj . 

The contributions proportional to cr in Eqs. (|36j) - (|41j) arise from the diagram in Figure 
1 (a) which has its origin from the U(1)a -D— term. We call this the trace contributions. For 
a non-anomalous U(l) gauge symmetry with universal scalar masses the trace term would 
vanish. However, for an anomalous U(l) gauge symmetry, trace of the flavor charges is 
not zero, so this term will induce flavor non-universal masses. The diagram in Figure 1 
(b) is the source of flavor non-universal contributions proportional to the gaugino mass 
M Af in Eqs. §^-g2>- 

Now we give the expressions for the one-loop contributions to the /3-function of the 
SUSY breaking A-terms of Eq. (0). Let us introduce the following notation: 

4 - «{,<'■'■■ (43) 
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Figure 1: (a) Trace correction from iD^-term to the soft scalar masses, (b) U(1)a gaugino-induced 
corrections to the soft masses. 



There are two types of contributions to the /3-functions of A^: one from the gaugino and 
the other from the A-terms. The flavor gaugino contribution arises from diagrams such 
as the one in Figure 121 The A-term contribution to (3 M/) cannot have the flavon field 
S propagating in the loop, so that contribution is included in the MSSM piece. 
The gaugino vertex contribution to (3 M| 3 -J (see Figure |2J) is 

^ a %) V = ^ M ^Vl (flf ^ + + # + n iM<li + + + \n%{n% - l)g s 2 ) . (44) 
Eq. (|44|) is obtained by summing all possible gaugino exchange diagrams. 




H. 



* S, n.. 



Fi gure 2: U(1)a gaugino-induced vertex correction diagram for A-terms. 



We now list the full one-loop (3 function for each AL. This generalizes the results of 
Martin EH- 



(3{A e 



f3(A< 



}MSSM 
>ij 



+ ^ {A e \f vX Y v - 2 ((gf f + {q]f + h 2 ) gl 
+ 2Y^A^ i3 + ±- 2 glZ^M XF , 



(45) 



16tt 



17T 2 i 



5Y ut Y u + Y et Y e + Tr f 3Y u Y ut + Y u Y ut 
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-3g 2 2 - 3^/5 - 2 ((gf ) 2 + (^) 2 + /> 2 ) 



-21^ 



+ W A e + Tr f 3A u F ut + A V Y V ^ 



+3M l9 */5 + 3M 2 <7 2 2 ] }. . + — g 2 F Z^M Xp , 



(46) 



(3(A u ) l3 = P(A^ SSM + T ^{A u [TT(Y^)-2(( q f) 2 + ( q J) 2 + h 2 )g 2 F 



+2Y«Tr {A^)} u + JL«£ Z«^M Af , 



(47) 



d\MSSM 



1 1 3 



^9 2 F A ij ((?? 



{q"f + V 



471-2 



Here we defined the combination of the U(1)a charges Z{, as 



q?q] + q?h + qjh + n% q s {q? + q] + h) + in£(n£ - l)q 2 s 



2 
1 



q?qj + q?h + qjh + n^qf + q d - + h) + ^<(n y 



From the charges listed in Table ^ we have 

/(ll, 13, 16) 
Z e = - (10,11,13) 
V (10, 11, 13) 



in Model 1 [a = 0) and 



((7,9, 12) 
(6, 7, 9) 
V (6, 7, 9) 



4,6,9) (1,3,6) 
3,4,6) (0,1,3) 
3,4,6) (0,1,3) 



4,6,9) (1,3,6) 
3,4,6) (0,1,3) 
3,4,6) (0,1,3) 



qtq) + q?h + q]h + n%q s (q^ + q] + h) + -n?-(n?- - l)q 2 s , 

q-q] + q*h + <$h + < i?s (^ + q] + h) + \k 3 {u\ 3 - i) q 2 . 



(4* 



(49) 



(50) 



(51) 



in Model 2 (a = 1) for the three different values of p = (0, 1, 2). 



4 Lepton Flavor Violating Decays /i — > &y and r — > /17 
4.1 Qualitative Analysis 

In the Standard Model with small neutrino mass, lepton flavor violating processes are 
highly suppressed. On the other hand, in the presence of low energy supersymmetry, 
LFV effects can be quite significant. In particular, LFV induced by the right-handed 
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neutrino Yukawa couplings in the MSSM can lead to p, — > c~) and r — > p'y decay rates 
near the current experimental limits 

Here we focus on flavor violation in leptonic processes. The slepton soft masses are 
more sensitive to the U(1)a gaugino corrections compared to those in the squark sector. 
This is because flavor violation in the squark sector is diluted by the gluino focusing effects. 
This is especially so when one considers the cosmological constraints on the lightest SUSY 
particle (LSP) mass. Demanding that the neutralino LSP constitutes an acceptable cold 
dark matter imposes the condition m — M1/2/4.4 in the context of supergravity models. 
This condition results from the coannihillation mechanism j22l for diluting the dark 
matter density which requires Tr mass to be about 5 — 15 GeV above the LSP mass. 

The approximate formulae for the sfermion soft masses in terms of the universal soft 
scalar mass mo and the common gaugino mass M\/i (for small to medium tan/3) are 

m 2 L ~m 2 + 0.52M 1/2 , 



m;~m z + 0.15M 1/2 , 

/2, 



m 2 n ~wiq + 6.5Mj 2 12 



m 2 u ~m 2 d ~m 2 + Q.lM 2 /2 . (52) 

From these expressions with m — Mi/ 2 /4.4 we see that the gaugino focusing effects make 
the squark soft masses much less sensitive to any flavor violating contributions. 

The right-handed neutrino induced LFV effects in our models depend on the overall 
factor e p in Eq. (JHJ). These processes will be suppressed for p — 1, 2 corresponding to low 
values of tan/3. As we noted in Sec. 2.3, p = 2 is preferred for leptogenesis, in which case 
the v c effects are small. 

There are three different sources of LFV in our models: (i) RGE effects between M st 
and the U(1)a symmetry breaking scale Mp induced by the U (1) gaugino, (ii) RGE effects 
between M st and the right-handed neutrino mass scale Mr induced by the neutrino Dirac 
Yukawa couplings, and (iii) the U(1)a D—term. Here we discuss only the RGE effects (i) 
and (ii). We call them the flavor gaugino induced LFV and z/ c -induced LFV. (The D-term 
contributions are not included in our numerical analysis, but are discussed in subsection 
4.3.) We give approximate formulas for these LFV processes by integrating the relevant 
/3-functions derived in Section 3. 

We adopt the minimal supergravity scenario (mSUGRA) for supersymmetry breaking. 
We assume universality of scalar masses and proportionality of the A-terms and the 
respective Yukawa couplings at the string scale. Gaugino mass unification is also assumed. 

The various flavor violating effects are summarized below: 

(1) Right-handed neutrino contributions to the scalar soft masses arising from Eq. 
(jHEJ) proportional to the Dirac neutrino Yukawa couplings: 

s{A l)^-(y^y % ^ + ^^^. (53, 

(2) Trace correction from D^-term in Eqs. (|36|) and (|37|) from Figure 1(a): 



5 (m 2 L ) A ~ -qf'lqslg 2 ^ 3m„ ^ (n£ + 

13 \ »=1,3 



xi - 2 
n m n — m„ 



i=l,3 



\n(M st /M F ) 
8vr 2 
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S (m 2 Y. ~ -qt\q s \g 2 F$ij 3 m 2 ^ ( n JJ + 

J v i=i > 3 



+™3 E K + O + ^ - *2 1 ln( ^{ MF) • (54) 



i=l,3 



(3) Gaugino mass correction from Figure 1(b): 

\ G ( l \ 2 r /w ^ 2 ln(M st /M F ) 



2tt 



5 (m 2 ) ~ <J y (M Af ) — . (55) 

(4) Right-handed neutrino induced vertex correction to the A e -terms (see Eq. (|45jl ): 

6At^-3A (Y^Y^) W^M , (56) 

(5) Flavor gaugino vertex correction to the j4 e -terms arising from Figure 2 (see the 
last term of Eq. (1351)): 



SAij ~ -m Af ^ f — . (57) 



In addition, we have flavor charge dependent wave function renormalization of the A— 
terms as given in Eq. 1)450. These are however not significant since they are diagonalized 
simultaneously with the corresponding Yukawa couplings. On the other hand, the vertex 
corrections to the A-terms given in Eqs. ((SH}) and (|57jl will induce nonproportionality in 
going from M st to Mp. 

The matrix elements Zfj in Eq. (|57j) are given in Eqs. (|5U|) and (|5T|) for different 
values of p. The elements in the (1,2) block of Z e are rather different from each other, 
suggesting that the gaugino vertex contributions can be very important for the process 
\i — ► &y. On the other hand, the elements in the second and the third rows are identical, 
hence, A\ z and A^ run at the same rate as their corresponding Yukawa couplings do in 
the short momentum interval. Therefore, this vertex correction for the process r — > fi^y is 
always suppressed in models with the texture of Eq. (jSJ). 

For fi — > e7 we find that the most dominant effect is from the flavor gaugino con- 
tributions to the soft masses. This is due to the following reason. It is proportional to 
the flavor charge squared and to the flavor gaugino mass squared (recall that we have 
m ~ Mx/2/4.4), both of which are large. On the other hand, the trace contributions to 
the soft masses depend linearly on the flavor charges and are proportional to uiq, which 
make them relatively small although the trace of the U(1)a charges itself is large. The 
right-handed neutrino contributions are significant only for p — 0. For other values of p 
the ^-contributions to the branching ratio for Zj — ► Ijj is suppressed by e 4p . 

We find that the gaugino contribution to the r — > /17 decay rate is always suppressed 
since tl and hl have the same flavor charges and since the tr-^r mixing angle is of order 
e 2 . The only significant effect to this process is from the right-handed neutrino effects 
when p = 0. 
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4.2 Numerical Results 



In this section we present our numerical results for the LFV processes p — > cy and r — > p*y. 
We adopt the mSUGRA scenario for the SUSY breaking parameters. At the string scale, 
taken to be M st = 10 17 GeV, we assume a universal scalar mass tuq and a common gaugino 
soft masses Mi/ 2 . The unified gauge coupling at 2 x 10 16 GeV is taken to be «g — 1/24. 
We assume the U(1)a gauge coupling g F to be equal to g 2 at the string scale. We evolve 
the soft SUSY breaking parameters from M st to the U(1)a gaugino mass Mp ~ M st /80 
(see Eqs. (|15|) and We use the numerical values of the Yukawa couplings given in 

Eqs. ip7) - Pty for this evolution. 

As explained previously, we take mo = M1/2/4.4 so that the relic abundance of neu- 
tralino dark matter can be reproduced correctly. With this choice we always find the 
neutralino to be the LSP with the tr mass higher than the LSP mass by 5 — 15 GeV. 
We impose radiative electroweak symmetry breaking condition. The SUSY higgs mass 
parameter p is chosen to be positive which is favored by b — > 57. 

We take M 1 / 2 to vary in the range 250 GeV to 1 TeV. The lower value satisfies 
the lightest higgs boson mass limit. We present the results for three different values 
of tan/3 = (5, 10, 20). The corresponding values of the exponent p are taken to be 
p = (2, 1, 0). The results are presented for two different values of A Q = (0, 300) GeV. 
When tan/? = 20, the lower limit on M1/2 is around 300 GeV, or else the radiative 
electroweak symmetry breaking would fail. The branching ratios are plotted against 
universal gaugino mass M1/2 in Figures 3-14. 

In Figure El we plot the branching ratio for the process p — > ey including all the LFV 
effects described earlier for Model 1 as a function of Myz for different values of tan/3 and 
A . We see that the branching ratio is in the experimentally interesting range for most of 
the parameter space. In this Figure for tan f3 = 20 we also show the branching ratio when 
the Dirac neutrino Yukawa coupling effects are maximized (denoted by "large Y v "). The 
horizontal line corresponds to the current experimental limit. 

In Figure |U the combined effect for p — > &y is plotted for Model 2. 

In Figure El we plot B (/i — > erj) induced solely by the right-handed neutrino Yukawa 
couplings. This result is identical for Models 1 and 2 since neutrino textures are the same 
for the two models. In Figure El we plot the branching ratio induced by the right-handed 
neutrino effects and the flavor gaugino effects for Model 1. Figure [7| has the same plot for 
Model 2. In Figure 8 (9) we plot B (p — > &y) induced by the trace term and the right- 
handed neutrino for Model 1 (2). Figure 10 (11) is a plot of the branching ratio including 
the effects of A-terms and v c for Model 1 (2). Figures IT21 and IT31 are the branching ratios 
for t — > ^7 including all LFV effects for Model 1 and Model 2. Figure IT41 which is valid 
for both Models 1 and 2, has the branching ratios for r — > pry induced only by the v c 
Yukawa coupling effects. 

From these figures we see that the decay p, — > &y is within the reach of forthcoming 
experiments. Discovery of r — >• p'y decay will strongly hint, within our framework, an 
origin related to the right-handed neutrino Yukawa couplings. 
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Figure 3: Branching ratio for the process fj, — * including all corrections for Model 1. The solid line 
corresponds to Aq = 300 GeV and the dashed line corresponds to Aq = GeV. For tan j3 = 20 we give 
two sets of curves, the upper one corresponds to the maximal value of the neutrino Yukawa coupling 
Y v . Here and in other plots, the straight horizontal line corresponds to the current experimental limit 
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Figure 4: Branching ratio for the process [i — > ej including all corrections for Model 2. 
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Figure 5: Branching ratio for the process — > induced by only the right-handed neutrino Yukawa 
coupling effects. This result holds for both Models 1 and 2. 
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Figure 6: Branching ratio for the process fj. — > induced by the gaugino corrections (plus v c effects) 
for Model 1. 
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Figure 7: Branching ratio for the process \i — > e~f induced by the gaugino corrections (plus v c effects) 
for Model 2. 
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Figure 9: Branching ratio for the process fj, — > e-f induced by the trace correction (plus v c effects) for 
Model 2. 
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Figure 10: Branching ratio for the process \i — > from the vertex corrections (plus v c effects) for 
Model 1. 
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Figure 11: Branching ratio for the process fi — ► e~f from the vertex corrections (plus v c effects) for 
Model 2. 



21 




400 600 800 1000 

M 1/2 (GeV) 



Figure 12: Branching ratio for the process r — » jjq including all the effects for Model 1. 




Figure 13: Branching ratio for the process r — * jjq including all the effects for Model 2. 
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Figure 14: Branching ratio for the process r — * /17 induced by only the right-handed neutrino Yukawa 
coupling effects for Models 1 and 2. 



4.3 The D— term Splitting Problem 

Any model based on a gauged flavor symmetry has a potential .D-term splitting problem, 
which could give rise to large FCNC processes even with a universal choice of soft scalar 
masses. Here we quantify this problem in the class of anomalous U(l) models. We 
point out that this problem is not as serious as it might naively appear, if the hierarchy 
mo = Mi/2/4.4 needed for an acceptable relic abundance of neutralino dark matter is 
assumed. 

The .D^-term contribution to the scalar potential including soft SUSY breaking mass 
for the flavon field S is given by 

V = ml\S\> + i«=M (JL. - + £ q^J (58) 

where q s and g« are the U(1)a flavor charges of the S field and the MSSM fields 
Minimizing the potential one finds the mass splitting among the MSSM sfermions to be 

~ 2 ~ 2 Qi ~ Qj ~ 2 /rn\ 

rrij. — rrif. = — : — p-m s . (59) 
I Qs I 

If a universal scalar mass is assumed even for this splitting could be unacceptably 
large. However, if we choose mo = Mi/2/4.4, this may be in an acceptable range. The 
low energy values for the slepton soft masses including the D^-term corrections to Eq. 
(1521) are 

(ml).. = r^ifal + m 2 + 0.52 A/? /2 ~ prm^ + 11.1 mjj, 
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(fhf).. = -^m 2 s + m 2 + 0.15 M 2 /2 ~ -^m 2 s + 3.9 m 2 . (60) 

Since the D-terms only contribute to the diagonal slepton mass splittings, any flavor 
violation arising from it must involve fermionic mixing angles. We find 

5ql 2 Y£ x m 2 s 



V /12 Y 2 e 2 (m 2 e ) 22 
( 5 L L y „ Sq^rh 2 

where {^ RR ^) 12 is the e,R-fi-R mixing parameter in the supersymmetric basis. Here 6qf 2 = 

(Qi ~ ?!)/|?*| = 2 - a and 5qf 2 = (gf - gfVI&l = !> with a = (1) for Model 1 (Model 
2). For — > e7 this gives 

(o; 2 - g« (2 ' a) r° ( ? 2 - < 63 > 

v J 12 1 j 1_ ( «k. ) 

1.95 V"X> / 

(^ LL ); - <? £ u ; a , (64) 

v y 12 -j j p_ / 2k. 1 

1 11.1 V»W 



where 



C. 



3.9 Y 2 \ V86.9' 46.9 
^i e 2 / 1 1 



C L = — ^ , — for (Model 1, Model 2). (65) 

II.IV22 V75.6 72. J v ; K ' 

From the experimental bound on \i — * cy one approximately has [8] 

For tan/3 = 5 and for Msusy ~ M1/2 (which is a reasonable choice) we find 

(S LL ' RR Y l2 \ exp < 0.6 x 10~ 3 (l.O x 10~ 2 ) for M V 2 = 250 (1000) GeV. (67) 
This gives the following constraint on m s jm^: 

Tfl 

— < 0.2 (0.9) for Mi/a = 250 (1000) GeV (68) 
m 

for Model 1 and 

— < 0.17 (0.7) for Mi/2 = 250 (1000) GeV (69) 
m 

for Model 2. For low values of Mi/ 2 this gives a significant constraint on rh s , while for 
large values of Mi/ 2 universality of all scalar masses may be maintained. 
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We wish to note that the m 2 appearing in Eq. (59) and in the subsequent discussions 
is the soft mass-squared of the scalar flavon field S evaluated at the U(1)a breaking 
scale Mp. Since there are singlet fields in the theory needed for anomaly cancellation 
which are active between the U(1)a scale and M st (Cf: Eq. (9)), the S field can have 
renormalizable Yukawa couplings of the type YkX^X^S where Xk are fields neutral under 
U(1)a- We have examined the evolution of m 2 between Mp and M st in the presence of 
such Yukawa couplings and found that even with the limits on m s given in Eqs. (68) and 
(69), m 2 = ml can be attained at the string scale for a range of Yukawa couplings Y k . 

The RGE running of m 2 between the string scale and the U(1)a symmetry breaking 
scale contains the same trace term and the same gaugino contributions that we considered 
for the MSSM sfermions in Eqs. (JHlj) and (JHHj) (see also terms proportional to o in Eqs. 
(36)-(41)) 



Sm 2 ~ -q s \q s \gp \3m 2 ]T (n£ + n% 

\ i=l,3 



J_ 2 V- / „x, X 2 ~ 2 \^(M st /M F ) 



i=l,3 



111 2 



,( ,2 (M ,2 In (M st /M F ) 

+ {q s 9F) {M Xf ) — . (70) 

2ir z 

Interestingly, this contribution from the .D-term evolution cancels the "trace contribu- 
tions" in the sfermion masses (the terms proportional to o in Eqs. (36)- (41)) |25j . 



5 Conclusion 

In this paper we have studied lepton flavor violation induced by a flavor-dependent anoma- 
lous U(l) gauge symmetry of string origin in a class of models which addresses the fermion 
mass hierarchy problem via the Froggatt-Nielsen mechanism. We have derived a general 
set of renormalization group equations for the evolution of soft SUSY breaking parameters 
in the presence of higher dimensional operators. These results should be applicable to a 
large class of fermion mass models. We have shown that the U(1)a sector induces signif- 
icant flavor violation in the SUSY breaking parameters during the RGE evolution from 
the string scale to the flavor symmetry breaking scale, even though this momentum range 
is very short. We have identified several sources of flavor violation: the U(1)a gaugino 
contribution to the scalar masses which is flavor dependent, a contribution proportional 
to the trace of U(1)a charge which is also flavor dependent, non-proportional A-terms 
arising from the U(1)a gaugino vertex correction diagrams, and the U(1)a -D-term. In 
addition, there are flavor violating effects in the charged lepton sector arising from the 
right-handed neutrino Yukawa couplings, which have also been included in our numerical 
analysis. The resulting flavor violation in the leptonic decays fi — >• e7 and r — > /ry are 
found to be in the experimentally interesting range. 

Adopting the minimal supergravity scenario for SUSY breaking, and choosing param- 
eters such that the needed relic abundance of neutralino dark matter is realized, we have 
presented results for the branching ratios B(fi — > ej) and B{r —>■ yrf) in two specific 
models of fermion masses. Figures 3 and 4 are our main results for the two models for 
B(/j, — > erf), while Figures 13 and 14 are our results for B(t — > y/y). The former should be 
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accessible to forthcoming experiments, while the latter is also in the observable range. Al- 
though we focused on two specific fermion mass textures these effects should be significant 
in a large class of models. 
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